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Abstract 


An analysis-based approach for developing shell-buckling design cri- 
teria for laminated-composite cylindrical shells that accurately accounts 
for the effects of initial geometric imperfections is presented. With this 
approach, measured initial geometric imperfection data from six graph- 
ite-epoxy shells are used to determine a manufacturing-process-specific 
imperfection signature for these shells. This imperfection signature is 
then used as input into nonlinear finite-element analyses. The imperfec- 
tion signature represents a "first-approximation" mean imperfection 
shape that is suitable for developing preliminary -design data. Compari- 
sons of test data and analytical results obtained by using several differen t 
imperfection shapes are presented for selected shells. Overall, the results 
indicate that the analysis-based approach presen ted for developing reli- 
able preliminary-design criteria has the potential to provide improved, 
less conservative buckling-load estimates, and to reduce the weight and 
cost of developing buckling-resistant shell structures. 


Introduction 

High-performance aerospace shell structures are inherently thin walled because of weight and perfor- 
mance considerations and are often subjected to destabilizing loads. In addition, fiber-reinforced compos- 
ite materials are becoming more common in aerospace vehicle structures because of their high strength, 
high stiffness, and amenability to structural tailoring that permits design solutions that are not possible 
with metals. Thus, buckling is an important and often critical consideration in the design of these struc- 
tures and reliable, validated design criteria for thin-walled shells are needed, especially for shells made of 
advanced composite materials. 
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As one might expect, shell-buckling design criteria have a history steeped in empiricism. ’ From 
approximately 1930 to 1967, many shell-buckling experiments were conducted on metallic shells. Typi- 
cally, the experiments yielded buckling loads that were substantially lower than the corresponding analyti- 
cal predictions, which were based on simplified linear bifurcation analyses of geometrically perfect shells 
with nominal dimensions and idealized support conditions. The seminal works by von Karman and 
Tsien 3 , by Donnell and Wan 4 , and by Koiter 5 identified small deviations from the idealized geometry of a 
shell, known as initial geometric imperfections, as a primary source of discrepancy between corresponding 
analytical predictions and experimental results. However, it appears that because of the intense computa- 
tional requirements of the nonlinear analyses used by von Karman and Tsien and by Donnell and Wan, and 
the asymptotic analyses used by Koiter, at that time, buckling design allowables were determined by estab- 
lishing lower bounds to test data. Specifically, empirical factors, that have become known as "knock- 
down" factors, were determined and were to be used in conjunction with linear bifurcation analyses for 
simply supported shells; that is, these empirical factors were used to "knockdown" the unconservative ana- 
lytical prediction. This approach to shell buckling design remains prominent in industry practice, as evi- 
denced by the extensive use of the NASA space vehicle design criteria. 6 ' 8 

From approximately 1960 to 1989, digital computers emerged and matured as a practical tool for sci- 
entists and engineers. With these computational tools, considerable effort was put forth to determine ana- 
lytically the effects of initial geometric imperfections on the buckling response of shells. These efforts 
were largely based on Koiter’ s asymptotic theory 1,5,9 for idealized structural shapes with nominal dimen- 
sions and also included investigations of shells made of fiber-reinforced composite materials. 10 ' 13 In addi- 
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tion, in-depth, substantial experimental investigations that focused on fiber-reinforced composite shells 
were conducted at institutions such as the NASA Langley Research Center and the Wright Aeronautical 
Laboratories at Wright-Patterson Air Force Base, in addition to various academic institutions worldwide. 
It was during this time period that the measurement of initial geometric imperfections in test specimens 
and high-fidelity response measurement techniques and apparatus were recognized as important elements 
in the development of design criteria for buckling of shells. 14-17 

Another significant development that occurred in parallel to the maturation of digital computers in the 
1960s to 1980s was the development of computational structural analysis codes based on finite-difference 
or finite-element methods. Two prominent examples that stimulated the development of this particular 
technology are the NASA Structural Analysis code (NASTRAN ) and the Structural Analysis of General 
Shells code (STAGS 19 ), which were funded by the United States government. The successors and inheri- 
tors of these codes now possess robust static and dynamic nonlinear analysis capabilities, along with 
extensive graphics capabilities, that enable scientists and engineers to model and visualize experiments or 
actual flight hardware to an enormous level of detail. Archetypal examples of the application of these 
tools are the in-depth, highly detailed studies of cylinder buckling response documented in Refs. 20-22 and 
the comprehensive and highly detailed nonlinear and buckling analyses of the Space Shuttle superlight- 
weight external fuel tank that is documented in Refs. 23-27. In addition, the extension of these analysis 
techniques and codes to other disciplines has helped to enable the development of advanced experimental- 
response measurement systems. Thus, the tools and understanding now exist for developing high-fidelity 
shell buckling design criteria for shells of general shape and material composition and for combined 
mechanical and thermal loads. 

Recent advances and recommendations for the development of modern shell-buckling design criteria 
have been given in Refs. 28 and 29. These references examine the potential benefits of adopting probabi- 
listic and hybrid probabilistic-deterministic approaches for developing design criteria. An important con- 
cept that is discussed in Ref. 29 and used in the present study is the establishment of characteristic 
imperfection signatures that are associated with the different shell-manufacturing processes. This concept 
is based on the idea that each manufacturing process will yield shells that possess essentially similar initial 
geometric imperfection shapes, with a well-defined range of pointwise amplitudes. With this information 
established, high-fidelity analyses and selective structural testing can be used to determine refined, reliable 
design criteria for shell buckling that are not overly conservative like the present lower-bound approach 
found in NASA SP 8007. 6 To this end, some effort has been put toward defining initial geometric imper- 
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fection models to be used in shell analysis for isotropic and composite shells. ’ Similarly, various prob- 
ability-based analysis methods have been proposed to develop reliable design buckling loads for imperfect 
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shells in which approximations to the measured imperfections are used in the analyses. ’ 

The present study focuses on the development of a refined shell-buckling preliminary-design criteria 
that is based on high-fidelity nonlinear finite-element analyses that include the effects of a manufacturing- 
process-specific geometric imperfection signature. Refined preliminary-design data are important 
because they substantially reduce the chance of getting "locked-in" to a poor design early in a vehicle- 
development program and, in the context of stability-critical thin-walled shell structures, can give a more 
accurate preliminary prediction of the load-carrying capacity of the shell. One main objective of the 
present study is to establish a "first-approximation" characteristic imperfection signature for selected cir- 
cular cylindrical shells made of a graphite-epoxy material that can be used to develop data suitable for pre- 
liminary design. Another main objective is to develop and validate an analysis-based shell-buckling 
design criteria that accurately accounts for the effects of initial geometric imperfections by using the 
imperfection signature defined herein. The final objective is to show the potential improvements in buck- 
ling predictions that can be realized, particularly for combined loading conditions. Towards these objec- 
tives, six graphite-epoxy shells with different shell-wall laminates that were previously tested and 
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analyzed are discussed first. Then, the nonlinear shell-analysis procedure that is used to predict the buck- 
ling and nonlinear response of the shells is described. Next, the details of how measured geometric imper- 
fections for six laminated graphite-epoxy shells are used to determine a manufacturing imperfection 
signature are presented. Then, analysis results and experimental results for selected compression-loaded 
shells are compared to validate the new high-fidelity analysis-based procedure for developing design crite- 
ria. Finally, numerically predicted buckling interaction curves are presented for three different shells sub- 
jected to combined axial compression and torsion to illustrate the potential benefits of the new procedure. 


Test Specimens, Imperfection Measurements, and Tests 

Six test specimens used in the present study are described in this section. These specimens were previ- 
ously fabricated and tested at the NASA Langley Research Center in support of a research program to 
study the effects of initial imperfections on the nonlinear response and buckling of compression-loaded 
composite shells. Much of this work has been documented in Refs. 20-22 and selected details about the 
test specimens, imperfection measurements, and tests are summarized in this section. Additional details 
on the cylinder tests can be found in these references. 

Test Specimens 

The specimens were fabricated from 12-in-wide, 0.005-in-thick AS4/3502 unidirectional graphite- 
epoxy tape material. The nominal unidirectional lamina properties of a typical 0.005 -in-thick ply with a 
fiber volume fraction of 0.62 are as follows: longitudinal modulus Ej = 18.5 Msi, transverse modulus E 2 = 
1.64 Msi, in-plane shear modulus G 12 = 0.87 Msi, and major Poisson’s ratio v ]2 = 0.30. The material was 
laid up by hand on a 15.75-in-diameter mandrel and vacuum bagged and cured in an autoclave to form six 
shells with different shell-wall laminates. These shells include shells with an 8-ply axially stiff [ +45 / 0 2 ] s 
laminate, an 8-ply circumferentially stiff [+45/90J s laminate, an 8-ply quasi-isotropic [ +45 / 0 / 90] s 
laminate, a 16-ply axially stiff [+45/0 2 ] 2s laminate, a 16-ply circumferentially stiff [+45/90 2 ], s lami- 
nate, and a 16-ply quasi- isotropic [+45 / 0 / 90] 2s laminate. These six shells are referred to herein as shells 
or specimens Cl through C6, respectively. These specimens had a nominal length of 16.0 in. and a nomi- 
nal radius of 8.0 in. The 8-ply and 16-ply specimens had a nominal shell-wall thickness of 0.04 in. and 
0.08 in., respectively, and had shell-radius-to-thickness ratios of 200 and 100, respectively. Both ends of 
the specimens were potted in an aluminum-filled epoxy resin to ensure that the ends of the specimen did 
not fail prematurely during the test. The potting material extended approximately 1.0 inch along the length 
of the specimens at each end, resulting in a test section that is approximately 14.0 in. long. The ends of the 
specimens were machined flat and parallel to facilitate proper load introduction during the tests. A photo- 
graph of a typical test specimen and the specimen coordinate system used to represent the corresponding 
geometry is shown in Fig. 1. The nominal cylinder length, test-section length, radius, and thickness are 
designated as L, L, R, and t, respectively. 

Imperfection Measurements 

Three-dimensional surveys of the inner shell-wall surface of the specimens were made prior to testing 
to determine the initial geometric shell-wall imperfection shape. Measurements were taken over a uniform 
grid with increments of 0.125 in. in the axial direction and 0.139 in. (approximately 1 deg of arc) in the cir- 
cumferential direction over the exposed surfaces of the specimens. A typical contour plot of the normal- 
ized initial geometric shell-wall mid-surface imperfections for specimen C3 is shown in Fig. 2. 

The measured shell-wall imperfection w 0 is normalized by the nominal shell-wall thickness t nom = 0.04 
inches. These results indicate that the initial geometric shell-wall imperfection is periodic in the circum- 
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ferential direction and exhibits slight variations in the axial direction. The amplitude of the imperfection 
varies from +1 ,21t nom to -1.536 t nom . 

Test Apparatus and Tests 

The specimens were instrumented with electrical-resistance strain gages and direct-current differential 
transducers (DCDT’s) were used to measure displacements. Three non-colinear DCDT’s were positioned 
in three corners of the upper loading platen and used to measure the end-shortening displacement and the 
platen rotations. A shadow moire interferometry technique was used to observe the shell-wall prebuck- 
ling, buckling and postbuckling radial deformation patterns. All data were recorded with a data acquisi- 
tion system, and the moire patterns were recorded photographically, on videotape, and with a high-speed 
digital video camera. 

The specimens were loaded in compression with a 120-Kip hydraulic universal-testing machine by 
applying an end-shortening displacement to the shell ends. To control the load introduction into the speci- 
men, the upper loading platen was aligned with the loading surface of the specimen before the test by 
adjusting leveling bolts in the corners of the upper loading platen until strains measured by selected strain 
gages on the specimens indicated a uniform axial strain distribution in the shell wall. The specimens were 
loaded until global collapse of the shells occurred. 


Finite-Element Models and Analyses 

The finite-element models and analysis methods used in the present study are based on experimentally 
validated models presented Refs. 21-22. In particular, the results presented in these references indicate 
that highly accurate buckling responses can be predicted for compression-loaded shells if initial geometric 
imperfections, material property variability, nonuniform thickness distribution, and other features are 
modelled to a high degree of accuracy. Such validated models that include the effects of these details are 
generally deemed necessary to develop the most accurate, analysis-based, high-fidelity design criteria. In 
line with the objectives of the current study, only the effects of the initial geometric imperfection are 
included in the numerical models. However, the model details, such as the mesh refinement and idealized 
boundary conditions, and the analysis methods are the same as those used in the previous studies and, as a 
result, the solutions are expected to be converged. In addition, the validated finite-element models used 
for compression-loaded shells are also used for the analysis of the torsion-loaded shells and the shells sub- 
jected to combined torsion and compression presented herein, and are expected to be converged because of 
the finite-element models extensive refinement. Additional details on the finite-element models and anal- 
ysis procedures are found in Refs. 21 and 22. 

Finite-Element Models 

The shells considered in this study were analyzed with the STAGS (STructural Analysis of General 
Shells) nonlinear shell-analysis code. 19 STAGS is a finite-element code developed for the static and 
dynamic analysis of general shells, and includes the effects of geometric and material nonlinearities in the 
analysis. A typical finite-element model of a specimen is illustrated in Fig. 1. The standard 410 quadrilat- 
eral element from the STAGS element library was used in the models which is a flat facet-type element 
that is based on the Kirchoff-Love shell hypothesis and the nonlinear Lagrangian strain tensor. The ele- 
ment nodes include three translational degrees of freedom and three rotational degrees of freedom. Large 
element rotations are accounted for by using a corotational algorithm. Typically, the elements of the 
finite-element mesh are approximately 0.2-in. by 0.2-in. square. However, some of the elements near the 


4 


boundary region of the shell are somewhat smaller in order to accurately represent the bending boundary- 
layer response near the end of the shell. Each element possesses four integration points, which are distrib- 
uted in such a way as to provide a modeling resolution of approximately 0.1 -in. by 0.1 -in square. This 
integration-point spacing is on the order of the measurement-point spacing used when measuring the initial 
geometric imperfections of the specimens. 

Compression-loaded shells and shells subjected to combined axial compression and torsion are ana- 
lyzed in the present investigation. To simulate the constraints provided by the potting material in a com- 
pression-loaded shell, the circumferential and radial displacements v and w were set equal to zero in edge 
regions of the shell as indicated in Fig. 1. Similarly, for shells subjected to combined axial compression 
and torsion, the radial displacement w was set equal to zero and the circumferential displacement v was 
applied to simulate a torsion load. 

Geometrically perfect and imperfect shells were analyzed in the present investigation. Nominal shell 
geometry, laminate thickness and lamina mechanical properties were used for the finite-element models 
for the geometrically perfect shells. The finite-element model of the geometrically perfect shells were 
modified to include the effects of various geometric imperfection shapes. In particular, linear-bifurcation- 
mode imperfections, measured shell-wall geometric imperfections, and simulated initial geometric imper- 
fections based on a "first approximation" manufacturing imperfection signature were considered. The var- 
ious initial geometric shell-wall mid-surface imperfection shapes were included in the finite-element 
models by introducing an initial out-of-plane perturbation to each node of the mesh. A typical finite-ele- 
ment model contained approximately 130,000 degrees of freedom. 

Nonlinear Analysis Procedure 

The STAGS finite-element analysis code uses both the modified and full Newton methods for its non- 
linear solution algorithms. The Riks pseudo arc-length path-following method 34 is used to continue a 
solution past the limit points of a nonlinear response. With this strategy, the incrementally applied loading 
parameter is replaced by an arc-length along the solution path, which is then used as the independent load- 
ing parameter. The arc-length increments are automatically adjusted by the program as a function of the 
solution behavior. The transient-analysis option in STAGS uses proportional structural damping and an 
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implicit numerical time-integration method developed by Park. Mass and stiffness damping factors are 
used in the proportional structural damping and are defined as a = 2nv and (3 = L) 2nv, respectively, 
where £ denotes the fraction of critical damping and is assumed to equal 0.15 for the shells considered 
herein. It is assumed that the dynamic buckling response of a shell is dominated by a particular frequency 
of vibration, v, and is defined as the critical frequency of vibration determined from a linear vibration anal- 
ysis in STAGS. These assumed damping parameters generally result in a slightly underdamped structural 
response during the transient analysis. Additional information on the transient analysis procedure can be 
found in Ref. 36. 

The prebuckling, buckling and postbuckling responses of the shells were determined by using the fol- 
lowing analysis procedure. The prebuckling responses were determined by using the geometrically non- 
linear quasi-static analysis capability in STAGS. The Riks pseudo arc-length path-following method in 
STAGS was used to compute the initial shell response until just before buckling occurred. The unstable 
buckling response of a shell was predicted by using the nonlinear, transient-analysis option of the code. 
The transient analysis was initiated from an unstable equilibrium state close to the limit point by incre- 
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menting the end displacement by a small amount. An initial time step of 1.0 seconds was used in the 
analysis, and the time step was automatically adjusted by the program as a function of the solution behav- 
ior. The transient analysis was continued until the kinetic energy in the shell had dissipated to a negligible 
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level, which indicated that the transient response had attenuated. Once the transient analysis had attenu- 
ated to a near-steady-state solution, the load-relaxation option of the code was used to establish the stable 
static equilibrium state. Results from conventional linear bifurcation buckling analyses were also deter- 
mined with STAGS for comparison with the nonlinear response results. 


Manufacturing Imperfection Signature 

Shell-wall geometric imperfections that had been measured previously (see Refs. 20-22) for the six 
laminated-composite shells considered herein were used to identify a manufacturing-process-specific "first 
approximation" imperfection signature for this type of structure. First, a reference or "best-fit" cylinder for 
the initial imperfection measurement data for each shell was determined. Specifically, a data reduction 
program that uses the method of least squares was used to compute the rigid-body translations and rota- 
tions and the mean radius of each cylinder. Then, the measured radial imperfections, deviations from a 
perfect cylinder, were recalculated with respect to the corresponding new best-fit cylinder. The measured 
imperfection data was then represented by a two-dimensional half-wave Fourier cosine series given by 



where R, L, and t nom are the shell radius, length and wall thickness; x and 8 are the axial and circumferen- 
tial coordinates; and k and / are integers denoting the number of axial half-waves and circumferential 
full-waves, respectively. Using a standard representation of this type for the measured imperfection 
enables convenient analysis and comparison of different imperfection distributions from different shells 
and different manufacturing processes. For example, a plot of the coefficient distribution of a half-wave 
Fourier cosine representation of the measured geometric imperfection shape for shell specimen C3 is pre- 
sented in Fig. 3. The results indicate that the largest-magnitude component of the imperfection is the k - 
0 and 1-2 component and that the magnitudes of the other components decrease rapidly as the wave 
numbers k and / increase. The k - 0 and 1-2 component corresponds to a cylinder with an elliptical 
cross-section that has a magnitude of 1 .2 t nom . The next largest component is associated with k = 1 and / 
= 2 and has a magnitude of approximately 0.2 t nom . In general, the dominant coefficients are associated 
with the long wave-length imperfections in the circumferential and axial directions. Coefficients for values 
of k and 1 >10 have a negligible effect on the overall imperfection shape. However, the short wave- 
length components are retained in the imperfection model since these modes may be associated with high 
imperfection sensitivity for the shells considered herein even though their magnitudes are generally very 
small. Results for shells Cl, C2, and C4-C6 indicate similar coefficient distributions and suggest that a 
manufacturing-process-specific imperfection signature for graphite-epoxy laminated-composite shells can 
be established in terms of a characteristic Fourier coefficient distribution and range. To this end, a mean 
imperfection shape and a corresponding standard deviation to the imperfection shape were determined 
based on the measured imperfection data. The mathematical expectation or mean value of wjx, 6) is 
given by 

q M .(A-,0) = (w o (x,0)) = X cos(^)| (A u ) cos (if-) + (B u ) sin(Z®e) J (2) 

where (...) denotes the mean value. A contour plot of the normalized mean geometric imperfections 
H w (x,6) is shown in Fig. 4. The mean shell-wall imperfection u H , is normalized by the nominal shell-wall 
thickness t nom = 0.04 inches. These results indicate that the initial geometric shell-wall imperfection is 
periodic in the circumferential direction, exhibits slight variations in the axial direction, and is similar in 
character to the measured imperfection shape for shell C3 shown in Fig. 2. The amplitude of the imperfec- 
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tion /u w varies from +1.77 t nom to -1.65 t nom . The corresponding standard deviation of the geometric 
imperfection shape <J w (x, 9) is defined by 

cr [w a C x,G )) = (( w 0 (x,6 ) - (w 0 (x,6))) 2 ) (3) 

A contour plot of the corresponding normalized standard deviation of the geometric imperfection shape is 
shown in Fig. 5. The results indicate that the maximum standard deviation is equal to 2.172 t nom and is 
primarily associated with long wave-length components of the imperfection shape. In addition, the coeffi- 
cient of variation was calculated for each component of the imperfection and gives the relationship 
between the mean and the standard deviation. The coefficient of variation, denoted by V, is given by 


V = 




where 

X k ,= \l< + B 2 kl 


(4) 

(5) 


A typical plot of V for an axial half-wave number k- 0 is shown in Fig. 6. The results indicate that V is 
relatively high for long wave-length imperfections; e.g., V - 1 .6 for k- 0 and 1-2, and decreases rapidly 
to a value of approximately 0.5 for / < 20. Results for the other values of k were also obtained and indi- 
cate similar trends. This information suggests that an approximate expression that relates the mean coeffi- 
cient values with the corresponding standard deviation can be developed. For example, an approximate V 
relationship for this particular case that is applicable for all values of k is given by 


V = 


1 . 25 - 2 ^/; 0</<20 
0 . 5 ; 20</<50 


(6) 


This equation is used in subsequent sections when defining "simulated" imperfections in place of the exact 
standard deviation function given in Eq. 3 to simplify the calculations. 


Response of Compression-Loaded Shells 

The effects of five different imperfection shapes on the response of the three compression-loaded 8- 
ply shells C1-C3 were investigated in the current study. Specifically, the effects of actual measured shell- 
wall imperfections, simulated imperfections that are based on the "first approximation" imperfection sig- 
nature, and linear-bifurcation-mode imperfections were analyzed. The simulated "first approximation" 
imperfection signature includes the mean imperfection with and without plus or minus one standard devia- 
tion. Load-shortening response curves and shell deformations are presented subsequently for the quasi- 
isotropic shell C3 that illustrate the overall response characteristics for the compression-loaded shell with 
the five different imperfection shapes. In addition, predicted buckling loads for the two orthotropic 8-ply 
shells C 1 and C2 with the five different imperfection shapes are also presented. Moreover, the predicted 
buckling loads are compared to the actual experimental buckling loads and to the design load one would 
get by extrapolating the results in NASA SP-8007 for isotropic shells to the composite shells. 

Four predicted load- shortening response curves are presented for shell C3 in Fig. 7a and include 
results for a shell with the actual measured imperfection shape; a shell with the mean imperfection shape, 
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denoted by |0. w ; a shell with the mean imperfection shape plus one standard deviation, denoted by p w + a w ; 
and a shell with the mean imperfection shape minus one standard deviation, denoted by p. w - G w . The 
axial load P and end-shortening A are normalized with respect to the linear bifurcation buckling load of the 
geometrically perfect nominal shell, P bif = 42,590 lbf, and the nominal shell-wall thickness, t nom = 0.04 in, 
respectively. The load-end-shortening curves clearly indicate a linear prebuckling response. The general 
instability occurs in the shell at normalized axial load levels that range from P/P bi f = 0.855, for a shell with 
p w + g w imperfection shape, to 0.943 for the shell with a p w - o w imperfection shape. The predicted 
buckling loads are summarized in Table 1. The general-instability response of the shell is followed by a 
sudden reduction in the axial load supported by the shell that is associated with the transient collapse 
response of the shell. The results indicate that the postbuckling load level is insensitive to the different 
types of imperfections considered herein. The corresponding load-time history of the transient collapse 
response of shell C3 with the actual measured imperfections are shown in Fig. 7b. The load-time history 
curve exhibits a sudden reduction in axial load until the collapse response attenuates and the axial load 
achieves a steady-state value. The kinetic energy in the shell obtains a maximum value during the tran- 
sient collapse response and dissipates over time and the shell reaches a stable postbuckling equilibrium 
state after approximately 0.007-0.008 seconds. The shell exhibits postbuckling load-carrying capacity, 
however, the effective axial stiffness of the specimen is significantly reduced in the postbuckling range of 
loading, as indicated by the large reduction in the slope of the load-shortening response curves. This 
reduction in effective axial stiffness is caused by large-magnitude radial deformations that develop in the 
specimen during buckling which result in significant load redistribution in the specimen and reduce the 
effective, load-carrying cross-section of the cylinder. 

The transient deformation responses for selected time steps during the transient collapse response of 
shell C3 with the actual measured imperfection, indicated by the letters A through D in Fig. 7b, are pre- 
sented in Fig. 8a through 8d, respectively. Just before buckling occurs, the shell wall deformations are 
characterized by several localized ellipse-like buckles, as indicated in Fig. 8a. The localization in the 
deformation pattern is caused by the combination of a local geometric shell-wall imperfection that is in the 
form of a significant variation in the shell-wall geometry and the bending boundary-layer deformations 
that form near the end of the shell. In addition, the localized deformations occur in regions with destabiliz- 
ing compressive axial and circumferential stresses. The coupling of these destabilizing stresses and the 
local deformations cause the buckling response in the shell. After approximately 0.0012 seconds have 
elapsed in the transient response, a single ellipse-like buckle has grown in amplitude and couples with the 
destabilizing stresses in the shell wall to cause the general instability and collapse of the shell. The magni- 
tude of the shell-wall radial displacement varies between +0.5 t nom . After additional time has elapsed in the 
transient collapse response, additional local buckles have formed around the circumference and along the 
length of the shell, as indicated in Fig. 8c, and the normalized axial load has decreased from 0.890 to 
0.554. In addition, some of the buckles in the shell begin to coalesce into larger diamond-shaped buckles. 
The magnitude of the shell- wall radial displacement varies between +3 t nom to -7 t nom . After approxi- 
mately 0.01 seconds have elapsed in the transient response, the kinetic energy in the shell has dissipated to 
a negligible level indicating that the transient response has attenuated, and the shell has deformed into a 
stable postbuckling mode-shape, as indicated in Fig. 8d. As loading continues in postbuckling region, the 
diamond-shaped buckles increase in size and the magnitude of the radial deformations of the buckles and 
the outer-surface ridges increase to between +4 t nom and -9 t nom . 

Deformation results were obtained for shell C3 with the mean imperfection shape, the p w + g w 
imperfection shape, and p w - o w imperfection shape and indicated similar response characteristics to 
those exhibited by the shell with the actual measured imperfection. In particular, predicted initial buckling 
deformations for shell C3 are presented in Figs. 9a through 9c. These results indicate that the initial buck- 
ling deformation response of a shell with the mean imperfection |i w and the | 0 . w + a w imperfection are 
similar to that exhibited by the shell with the actual measured imperfection and are characterized by a 


localization of a single buckle near one end of the shell. In contrast, the shell with the |0. w - o w imperfec- 
tion shape exhibits distributed initial buckling deformations shown in Fig 9c that are similar to those 
exhibited by the corresponding geometrically perfect shell. These results are significant because they sug- 
gest that the use of a mean imperfection shape and the corresponding standard deviation can bound the 
experimental buckling loads as well as predict physically meaningful response characteristics such as the 
buckling deformations. Results were also obtained for the two 8-ply orthotropic shells Cl and C2. In par- 
ticular, buckling loads for the various imperfection shapes were calculated and are also presented in Table 
1. Most notable, is the fact that the predicted buckling loads for the shells with the actual measured imper- 
fection are bounded by the results for the shells with the (i w + ct w imperfection shape, and ji w - G w 
imperfection shape. In addition, the experimental results agree well with the predicted results that include 
the effects of the actual measured geometric imperfection. The correlation would improve still further if 
other measured imperfections were included in the finite-element model; e.g., measured thickness varia- 
tion, material property variation, and loading nonuniformities, as shown in Refs. 21 and 22. Specifically, 
the measured thickness imperfection has been shown to account for an additional 5-10% reduction in the 
buckling load. In contrast, the results for shells with linear-bifurcation-mode imperfections are, for an 
imperfection amplitude equal to 0.5 times the nominal wall thickness, ultraconservative. Similarly, the 
lower-bound design recommendation given in NASA SP-8007 for the corresponding isotropic shell is 
equal to P cr / P^f = 0.47, which is also ultraconservative, as expected. 


Response of Shells Subjected to Combined Axial Compression and Torsion 

A numerical study of the effects of the four different imperfection shapes on the buckling load of the 
8-ply shells C1-C3 subjected to combined axial compression and torsion was also conducted. Normalized 
buckling interaction curves that include curves for a geometrically perfect shell, a shell with the actual 
measured imperfection, a shell with a q w + c w imperfection, a shell with a p w - o w imperfection, a shell 
with an linear-bifurcation-mode imperfection, and the corresponding NASA SP-8007 design curves are 
shown for shells C1-C3 in Figs. 10-12, respectively. The measured experimental buckling load for the 
compression loading is indicated on the figure by the solid circular symbol and is included for comparison. 
The buckling loads (stress resultants) N x and N xy are predicted from a geometrically nonlinear analysis, 
and are normalized with respect to the corresponding the linear bifurcation buckling loads N x bif , and 
N xy bi f of a geometrically perfect shell. The superscript "+" indicates results for a positive torsion load. A 
summary of linear bifurcation buckling loads for shells C 1 through C3 is presented in Table 2. 

The predicted results shown in Figs. 10 through 12 indicate that the buckling interaction curves for the 
geometrically perfect and imperfect shells are skewed with respect to the plane given by N xy — 0. The 
skewing is caused by the fact that the diagonal compression caused by the shear load is aligned with the 
outer 45-deg fibers for the case of negative shear and is perpendicular to the fibers for the case of positive 
shear. Thus, for a negative-shear load the shell is stiffer. In addition, the buckling interaction curves are 
concave and are highly dependent on the laminate stacking sequence. This behavior is similar to the buck- 
ling behavior of anisotropic plates that are subjected to combined compression and shear loads presented 
in Refs. 38 and 39. In contrast, the NASA SP-8007 design guideline assumes a linear relationship between 
the critical compression-buckling load and the torsion-buckling load. Moreover, the bilinear curve is 
symmetric about N xy = 0 because the results in NASA SP-8007 are for isotropic shells. 

The results in Figs. 10 through 12 show that shells Cl through C3 are sensitive to the as-measured 
imperfections and the p w + c? w imperfection when the shell is subjected to compression-dominated load- 
ing. However, as the torsion load increases, the predicted buckling loads exhibit significantly less imper- 
fection sensitivity for these particular imperfection shapes. The regions of the buckling load interaction 
curves for which the imperfection sensitivity exists is strongly dependent on the laminate stacking 
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sequence. For example, the buckling load interaction curves for shell Cl exhibit a relatively large region 
of imperfection sensitivity for normalized torsion loads N xy /N^ bif that ranges from -0.9 to 0.5 and nor- 
malized compression loads N x /N xbif that range from 0.65 to 0.97, as shown in Fig. 10. However, the 
regions of imperfection sensitivity in shells C2 and C3 are significantly smaller than that exhibited by shell 
Cl, as shown in Figs. 11 and 12, respectively. In contrast, the results indicate that a shell with a linear- 
bifurcation-mode imperfection exhibits a significant reduction in the buckling loads for all load cases. The 
large difference in the response trends is caused by a fundamental difference in the shape of the imperfec- 
tions in the shell and how they couple with the initial buckling deformation response of the shell. Specifi- 
cally, the linear-bifurcation-mode imperfection is similar in character to the initial buckling response in the 
shell and, thus, exhibits a very strong coupling with the initial buckling deformation of the shell. In con- 
trast, the actual measured imperfection and the imperfection signature, with imperfection amplitudes as 
large as w ( /t = 1.92, are primarily a long-wave imperfection and these imperfection shapes do not exhibit 
an artificial strong coupling with the buckling deformations in this shell; thus, the reduction in the buck- 
ling load imperfection sensitivity is much less pronounced. In general, the results for shells with linear- 
bifurcation-mode imperfections are, for the relatively small value of the imperfection amplitude consid- 
ered here, ultraconservative. Similarly, the lower-bound design recommendation obtained by using NASA 
SP-8007 for the corresponding isotropic shell is also ultraconservative, as expected, and misrepresent the 
behavioral trends. 

Finally, it is important to point out that although no test data were obtained for a C3-type shell sub- 
jected to torsion or combined torsion and compression, it is expected that the results for the shell with the 
actual measured imperfection will be very accurate. This expectation is based on the fact that torsion- 
loaded shells are not compression dominated like the corresponding compression-loaded shells and, as a 
result, will not exhibit as high of a degree of imperfection sensitivity. 


Design Implications 

The approach presented herein for developing shell-buckling design criteria appears to offer a means 
for obtaining preliminary-design estimates of buckling loads for laminated-composite cylindrical shells 
subjected to combined loads, without having to resort to extensive testing programs. In addition, this 
approach offers a relatively affordable alternative to relying on historical test data for shells that do not 
represent the configuration, material system, or fabrication process associated with the design of interest. 
This design-criteria-development approach could be used as a tool in the early stages of a design process to 
determine refined estimates of the buckling loads for the actual as-built shell structure under consideration. 
To insure success for problems of interest where refined design criteria are nonexistent, this approach 
should always be used with a selected number of carefully conducted tests that should be used to validate 
the finite-element analyses and provide confidence in the design results. The validation tests should be 
carefully conducted so that representative analysis results can be generated to correlate with test results. 
The provenience and pedigree of each test specimen should be determined or measured to assure that the 
appropriate shell-wall geometry and thickness distribution, material properties, and fabrication process- 
specific effects are included in the analysis used to develop the design criteria. 

The design-criteria-development approach described herein has wider implications; that is, it can be 
used to generate buckling design curves for commonly used structural forms and material systems that 
account for manufacturing-process imperfections and provide reliable estimates of buckling resistance that 
are suitable for preliminary design. Design curves of this fidelity will essentially provide a scientifically- 
based replacement for the overly conservative empirical knockdown factors typically used in design today 
to account for imperfection sensitivity of shell structures. Once a database of manufacturing-process-spe- 
cific imperfections is established by the technical community, a designer will be able to determine which 
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various imperfections are representative of the manufacturing procedure or fabrication process that is 
appropriate for his or her design and then use the relevant design curves. This approach will enable indus- 
try to develop a given design without having to repeatedly pay for extensive development tests to develop 
design knockdown factors. In addition, the refined design criteria should yield substantial weight savings 
in a given design. 


Concluding Remarks 


A new analysis-based approach for developing shell-buckling design criteria for laminated-composite 
cylindrical shells made of graphite-epoxy material with a tape-laying process has been presented. This 
approach accounts for the effects of the actual measured initial geometric imperfections in a manner that is 
suitable for the development of preliminary design data for this class of shells. In particular, a characteris- 
tic manufacturing imperfection signature has been identified from the actual measured geometric imper- 
fections for six laminated graphite-epoxy shells and the process used to synthesize the imperfection 
signature has been discussed in detail. More specifically, a "first-approximation" mean imperfection shape 
and the corresponding standard deviation of the imperfection shape were determined from the measured 
imperfection data for the six shells. Then, a "first-approximation" imperfection signature was used with 
high-fidelity nonlinear finite-element analyses to determine predictions of the buckling behavior of the 
three shells for uniaxial compression loading, torsion loading, and for combined uniaxial compression and 
torsion loading. Detailed comparisons of the analytical buckling predictions obtained by using the "first- 
approximation" imperfection signature, by using the actual measured imperfection data, and by using a lin- 
ear-bifurcation-mode imperfection shape, and with test results have been presented. Moreover, compari- 
sons with the corresponding results obtained by extrapolating the lower-bound design data given in NASA 
SP 8007 to the laminated-composite shells have been presented. In addition, the high-fidelity nonlinear 
shell analysis procedure used to predict the buckling and nonlinear response of the shells has been 
described and response characteristics have been discussed. 


The results presented herein indicate that the predicted buckling loads for the shells obtained by using 
the actual measured imperfection shapes in the finite-element analyses agree well the corresponding exper- 
imentally measured buckling loads for the compression-loaded shells, e.g., 10% maximum difference 
between predicted buckling loads with actual measured imperfection and the corresponding measured 
buckling loads. With regards to the "first-approximation" imperfection signature, the results show that, 
for the most part, the predicted buckling interaction curves for combined compression-torsion loadings 
that were obtained by using the actual measured imperfection shape fall within the corresponding interac- 
tion curves for plus or minus one standard deviation. In addition, comparisons of the nonlinear finite-ele- 
ment analyses show that the "first-approximation" imperfection signature yields response characteristics 
that are similar to those obtained by using the actual measured imperfection data in the analyses. These 
results are significant because one-standard-deviation bounds to the test data and the reproduction of sig- 
nificant nonlinear behavioral characteristics represent necessary conditions that should be met by a "first- 
approximation" imperfection signature that is to be used to develop an analysis-based design criteria. In 
contrast, the results also show that buckling interaction curves obtained analytically by using a linear- 
bifurcation-mode imperfection shape and the curves obtained by using the lower-bound experimental data 
in NASA SP 8007 are ultraconservative, compared to the corresponding results obtained analytically by 
using the actual measured imperfection data or the "first-approximation" imperfection signature. Overall, 
the results indicate that the analysis-based approach for developing reliable preliminary-design criteria that 
has been presented herein has the potential to provide improved, less conservative buckling-load estimates, 
and to reduce the weight and cost of developing buckling-resistant shell structures. 
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Table 1: Numerically Predicted and Experimentally Measured Buckling Loads for Compression-loaded Shells. 





Predicted Buckling Loads, lbf 




Shell 

Perfect 

Shell 

Measured 

Imperfection 

Mean 

Imperfection, 

R 

p + 0 

Imperfection 

p - 0 

Imperfection 

Linear- 

bifurcation- 

mode 

Imperfection, 
w <Anom ~ O--' 

Experimental 

Buckling 

Load 

Cl 

29900 

29300 

29000 

27500 

29850 

17040 

27788 

C2 

36300 

34400 

36170 

34430 

36930 

14470 

31924 

C3 

41500 

38180 

39100 

36400 

40200 

21300 

34082 


Table 2: Numerically Predicted Linear Bifurcation Buckling Loads (Stress Resultants) for Geometrically Perfect 8- 


ply Shells. 

Shell 

N xMP lbf/in 

N x ,bif’ lbf/in 

Cl 

592.4 

205.0 

C2 

729.4 

186.3 

C3 

847.2 

193.4 
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z, w 



Fig. 1 Typical specimen and finite-element model geometry and boundary conditions. 


w 0 =w 0 /t 


nom 
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1.272 

1.072 

0.871 

0.671 

0.470 

0.269 

0.069 

- 0.132 

- 0.333 

- 0.533 

- 0.734 

- 0.934 

- 1.135 

- 1.336 

- 1.536 


Circumferential location 0, degrees 


Fig. 2 Typical measured inner-surface imperfection shape for shell specimen C3. 
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Fig. 3 Coefficient distribution of a Fourier series representation of the measured imperfection shape for 

shell specimen C3. 


- 0.5 


- 0.25 


Axial 

location, 0 
x/L t 


0.25 


0.5 


n w =njt 


nom 



I 


1.773 

1.528 

1.283 

1.039 

0.794 

0.550 

0.305 

0.060 

- 0.184 

- 0.429 

- 0.673 

- 0.918 

- 1.162 

- 1.407 

- 1.652 


0 90 180 270 360 

Circumferential location 9, degrees 


Fig. 4 Calculated "first-approximation" mean imperfection signature. 
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® t nom 

2.172 
2.029 
1.886 
1.743 
1.600 
1.457 
1.314 
1.171 
1.027 
0.884 
0.741 
0.598 
0.455 
0.312 
0.169 

0 90 180 270 360 

Circumferential location 9, degrees 



Fig. 5 Calculated standard deviation to the imperfection signature shown in Fig. 4. 



number, k= 0 
variation 


Fig. 6 Coefficient of variation for the axial half-wave number k = 0 as a function of circumferential wave 
number /. 
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End-shortening, M nom 



a) Predicted load-end-shortening response curves for b) Predicted load-time history during collapse for 
shell C3 with selected imperfections. shell C3 with the actual measured imperfection. 


Fig. 7 Predicted nonlinear response for the compression-loaded quasi-isotropic shell C3 with various 
imperfection shapes. 



a) Prebuckling deformations b) Initial buckling 
Time = 0.0 seconds deformations 

Time = 0.0012 seconds 
P/P„,. = 0.890 


P/P cr = 0.896 


c) Transient buckling 

deformations 

Time = 0.0024 seconds 

P/P rr = 0.554 


d) Initial post-buckling 
deformations 
Time = 0.01 seconds 
P/P cr = 0.311 


Fig. 8 Predicted buckling deformation response for shell C3 with the actual measured initial imperfection. 
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a) Mean imperfection, p w 
Time = 0.001 1 seconds 


b) |i w + o w imperfection 
Time = 0.0012 seconds 


c) p w - O w imperfection 
Time = 0.0010 seconds 


P/P cr = 0.918 


P/P cr = 0.855 


P/P cr = 0.943 


Fig. 9 Typical initial buckling deformation responses for shell C3 with selected simulated initial imperfec- 
tions. 



AT 

,w xy,b/f 


Fig. 10 Predicted buckling interaction curves for the axially stiff orthotropic shell Cl subjected to com- 
bined axial compression and torsion. 
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,w xy,b/f 


Fig. 1 1 Predicted buckling interaction curves for the circumferentially stiff orthotropic shell C2 subjected 
to combined axial compression and torsion. 



,w xyb/f 


Fig. 12 Predicted buckling interaction curves for the quasi-isotropic shell C3 subjected to combined axial 
compression and torsion. 
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